Triple-g octupolar ordering in Np02 
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We report the results of resonant X-ray scattering experiments performed at the Np M^s edges 
in Np02. Below To = 25 K, the development of long-range order of Np electric quadrupoles is 
revealed by the growth of superlattice Bragg peaks. The electronic transition is not accompanied 
by any measurable crystallographic distortion, either internal or external, so the symmetry of the 
system remains cubic. The polarization and azimuthal dependence of the intensity of the resonant 
peaks is well reproduced assuming anisotropic tensor susceptibility (ATS) scattering from a triple-g 
longitudinal antiferroquadrupolar structure. Electric quadrupole order in NpC>2 could be driven by 
the ordering at To of magnetic octupoles of Ts symmetry, splitting the Np ground state quartet and 
leading to a singlet ground state with zero dipole magnetic moment. 



PACS numbers: 75.25.+z,75.10.-b 

For half a century the low temperature properties of 
Np02 have mystified theorists and experimentalists alike. 
Upon cooling from room temperature a single phase tran- 
sition is observed at To ~ 25.5 K Similarities with 
UO2 suggested a magnetic nature of the phase transi- 
tion, but no magnetic order was found by neutron diffrac- 
tion [|| nor by MoBbauer spectroscopy Q, which estab- 
lished an upper limit of « 0.01 /l*b for the ordered mag- 
netic moment /2 rd- However, Np + ions in NpC>2 are 
Kramers ions (5/ 3 , 4 19/2)- In the absence of interac- 
tions breaking time-reversal symmetry the ground state 
has to carry a magnetic moment, pi, whatever the crys- 
talline environment. A fluctuating magnetic moment of 
finite size would be revealed by a Curie-like divergence 
of the susceptibility at low temperatures, whereas the 
experiments reveal a flat susceptibility between 15 and 
5 K [||. Moreover, no evidence for a crystallographic dis- 
tortion, neither external nor internal, has been found by 
synchrotron experiments Q . 

Another clement of the puzzle was recently provided by 
muon spin relaxation experiments, showing the abrupt 
appearance of a precession signal below T 0. This 
implies that the order parameter (OP) sets up a mag- 
netic field at the muon stopping site and provides defini- 
tive evidence that the OP breaks invariance under time 
reversal. By assuming antiferromagnetic (AF) order of 
the same kind of that established for UO2, i.e. a type-I, 
triple-(f structure, the authors deduced an ordered mo- 
ment fi OT d ~ 0.1 /ie, a value much larger than the upper 
limit compatible with MoBbauer spectroscopy. In par- 
allel with this finding, direct evidence for long-range or- 
der in Np02 was obtained through the observation below 



To of superlattice reflections in resonant X-ray scattering 
(RXS) experiments at the Np M^s absorption edges [||. 
The superstructure Bragg peaks occur at Q = G+ (0 1) 
positions, where G is a reciprocal lattice vector. This is 
the same periodicity found for the AF phase in UO2, 
and the observations were taken as evidence for the oc- 
currence of longitudinal triple-<f AF order. Santini and 
Amoretti ||] pointed out the possibility of explaining the 
whole body of experimental evidence assuming magnetic- 
octupole order instead of magnetic-dipole order. Octupo- 
lar order would lift the degeneracy of the r§ Np ground 
state and generate an interstitial magnetic field, in agree- 
ment with neutron spectroscopy ]^| and muon spin reso- 
nance results. However, octupolar order can be directly 
observed in RXS only through E2-resonances, whilst the 
resonances observed in Q are at the El absorption edge. 
Indeed, all previous experimental data taken at the ac- 
tinide M^s edges indicates a very strong dominance of 
El-processes (3d 3 / 2 ,5/2 «-> 5/). To our knowledge there 
is no evidence for E2 contributions (3d 3 / 2 ,5/2 6*7), 
and indeed matrix elements involving E2 promotion to 
63 states are expected to be too small to allow for any 
observable signal. 

To clarify the above confusion, we have undertaken a 
new RXS experiment at the magnetic scattering beam- 
line, ID20, of the ESRF. We performed polarization anal- 
ysis of the diffracted radiation and measured the depen- 
dence of the intensity from the azimuthal angle, ifi (the 
angle describing the rotation of the crystal about the 
scattering vector). These analyzes were not possible in 
the experiment reported earlier in || , preventing an un- 
ambiguous determination of the origin of the resonance. 
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FIG. 1: Energy scans through the Np M4 absorption edge. 
The o — > 7r (main panel) and a — > a (not shown) resonances 
of the (0 3) superstructure peak have a Lorentzian-squared 
(solid line) rather than a Lorentzian (dashed line) line shape. 
The maximum at the resonance lies at the absorption edge as 
measured at the (0 2) Bragg reflection (insert). 



Indeed, the results we present in this Letter show that 
the superlattice peaks in Np02 are not magnetic, but 
arise from the asphericity of the Np 5/ electron density 
leading to an anomalous tensor component in the atomic 
scattering factor. In other words, the superlattice peaks 
signal the occurrence of electric-quadrupole long range 
order below To. This conclusion is incompatible with 
the particular octupolar model given in Ref. ||, which 
predicts an undistorted charge density for the 5/ ground 
state, with vanishing quadrupole moment. 

The experiment was performed on a single crystal of 
0.7 x 0.7 x 0.2mm 3 in volume, with a flat (0 1) surface. 
A closed-cycle refrigerator equipped with an azimuthal- 
rotation stage provided a base temperature of 12 K. An 
Au(l 1 1) single crystal was used to analyze whether the 
polarization of the scattered beam was parallel (n) or 
perpendicular (a) to the scattering plane. The incident 
beam had a polarization and the scattering plane was 
vertical. 

Fig. |l] shows the a — ► 7r intensity of the Q = (0 3) su- 
perlattice reflection as a function of the photon energy, 
E, around the Np M4 absorption edge. The data can 
be fit to a Lorentzian squared line-shape, centered near 
the El threshold. The presence of a resonance in the 
o~ — > <j channel (not shown) excludes that the scattering 
arises from dipole-magnetic order, as this would give only 
o~ — ► 7r resonant scattering. Studies measuring magnetic 
order at actinide M edges have shown Lorentzian energy 
profiles, whereas we clearly observe a Lorentzian-squared 
form, as it is expected when the intermediate state split- 
ting is much smaller than the core hole life time jlo| . 

The results of -0 (azimuthal) scans about the (0 3) 
peak, taken at T = 12 K with E = 3.846 keV, are shown 
in Fig. |^(a) for both the <r — > c and the a — > 7r channel. 
Fig. ||(b) and ^(c) show, respectively, the ifj dependence of 
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FIG. 2: Dependence of the scattering on the azimuthal angle, 
tp: (a) Intensity at the (0 3) quadrupolar reflection in the 
a- (filled circles) and it- (open circles) channels with incident 
polarization a. Different scaling factors have been applied to 
correct for differences in the efficiency of the PA in the a and n 
settings, (b) Stokes parameter Pi for the (0 1) (open trian- 
gles) and (0 3) (filled triangles) quadrupolar reflections, (c) 
Stokes parameter P2 for the (0 3) peak. The lines represent 
model calculations based on eqs. Q, and ^ 

the Stokes parameter Pi = (I a -^ a —I a ^ 1T )/(I a ^ a -+I cr ^ 1T ) 
and P 2 = (1+ - 0, for the (0 1) and 

(0 3) reflections. P2 can be considered as a measure 
of the phase relation between the two polarization chan- 
nels and is defined by the intensities measured with the 
analyzer oriented at ±45° with respect to the scattering 
plane, 

7± cx ±=\F a ^ a ± F a ^\ 2 (1) 

where F a ^ a and F a ^^ are the o — * a and a — > tt scat- 
tering amplitudes. 

The dependence of polarization and intensity of the 
scattered radiation on ip can be modeled with satisfac- 
tory agreement by assuming a triple-<f antiferroquadrupo- 
lar ordering of (T5) quadrupoles, and using the ATS cross 



section [[12| |lj] as detailed below and shown by the 
lines drawn in Fig. ^. 

Above To, NpC>2 crystallizes in the Fluorite structure 
with space group (SG) Fm3m, Np 4+ at 4a and O 2 " ions 
at 8c Wyckoff positions. A periodic distribution of elec- 
tric quadrupoles can be Fourier expanded, whatever the 
type of order. A triple-cf structure is obtained when 3 
components of the star of the propagation vector q en- 
ter in the Fourier sum. For <fi = (1 0), q\ = (0 1 0) 
and gjj = (0 1), the arrangement schematically shown 
in Fig. ^ is obtained, with charge distribution distorted 
along the (1 1 1) directions of the cubic unit cell. The 
SG symmetry is lowered from Fm3m to Pn3m, the only 
maximal non-isomorphic subgroup of Fm3m that is non- 
symmorphic and simple cubic. Within this SG, Np ions 
can be accommodated on the 46 positions, with reduced 
point symmetry D^d but the same crystallographic ex- 
tinction rules of the FCC para-quadrupolar SG. Oxy- 
gen ions can occupy two inequivalent positions (2a and 
6d), where all coordinates are uniquely fixed by symme- 
try alone, so that this electronic phase transition does 
not allow a shift of the oxygen ions. The symmetry of 
the ordered state remains cubic and, apart from a pos- 
sible change of the lattice parameter , the transition 
will not be accompanied by any distortion. In that case, 
the quadrupolar OP cannot be measured with neutron 
or conventional X-ray diffraction techniques. 

RXS from quadrupolar order has recently been ob- 
served in several systems |l6). It is well described 
within the framework of ATS scattering [0| [b|, 
which occurs when the photon energy is tuned to an ab- 
sorption edge of an atom. The anisotropy of this atom's 
polarizability may lead to a finite scattering cross section 
at reflections which are normally forbidden due to glide 
plane or screw axis extinction rules. The scattering am- 
plitude arising from El transitions can be described by 
second rank tensors, which are invariant under the point 
symmetry of the scattering atom. 

The scattering amplitude for the induced quadrupolar 
order must be represented by a symmetric tensor f(Q). 
For Q = (0 L), with L odd, the scattering length is 
given by F(Q) = t ■ f(Q) ■ e = $(e' x e y + e' y e x ), where 
e and e* are the polarization vectors of the incident and 
scattered beam, and $ is proportional to the Fourier com- 
ponent of the quadrupolar operator, <i>. With incident a 
polarization and ho = (1 0) as azimuthal reference vec- 
tor, we find 



F a ^ a = $sin(2V>) 
F a ^ = $sin(0) cos(2V>), 



(2) 
(3) 



where is the Bragg angle. The resulting scattered 
intensities are I a -+ a oc IF^o-I 2 and I a ^-K oc li 7 ^^ 2 . 
As shown in Fig. ||(a), the experimental data are well 
matched by this model with one overall scale factor. This 
factor can be eliminated by calculating the Stokes pa- 




FIG. 3: Crystal structure of Np02 in the antiferroquadmpo- 
lar state with space group Pn3m. The ellipsoids represent 
the orientation of the local symmetry axis at the Np position, 
not the actual charge distributions. The O atoms are shown 
as spheres. 



rameter P\. The data are shown in Fig. |^(b), along with 
model calculations. 

The quantities I± defined in eq. [l] are given by 



I± oc 



(1 - cos 2 (0) cos 2 (2i/;) ± sin(0) sin(4^)) . (4) 



The Stokes parameter P2 can be calculated, and is in 
good agreement with the experiments, Fig. ^(c). 

However, the quadrupolar order alone is not a suffi- 
cient ingredient, as it cannot explain the absence of a dis- 
ordered magnetic moment B and the breaking of invari- 
ance under time reversal . The lowest-rank multipolar 
OP consistent with the experimental findings (no ordered 
and fluctuating dipoles, broken time reversal symmetry) 
is an octupole P] . Under octahedral symmetry, the seven 
octupolar operators belong to irreducible representations 
T2, T4 or T^. A r4-octupole OP must be ruled out. 
In fact, this would always be accompanied by an order- 
ing of dipoles as the latter belong to T4 as well. A T2 
OP is consistent with most properties of Np02 ||, but 
it cannot explain the RXS results, as it does not carry 
an electric quadrupole moment and the a — * a signal 
would vanish Jl^j . T5 octupolar operators are given by 
symmetrized combinations of Oi = Ji(Jj — J|), where 
ijk = xyz, yzx or zxy. The little cogroup of the order- 
ing wavevector is D^h ■ ^5 decomposes into the tetragonal 
representations I4 (1-d) and (2-d), which identify 
two possible type-I octupolar orders, a "longitudinal" one 
and a "transverse" one. For a single-g structure, both 
types of OPs are inconsistent with the observed quench- 
ing of dipoles and with the present RXS results. There 
is just one single type of order which can explain all ob- 
servations, and this is a triplc-cf longitudinal structure, in 
which the OPs associated with the three wavevectors 
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of the star of q have the same amplitude, p. The simplest 
conceivable mean-field (MF) Hamiltonian for each of the 
four sublattices s includes the crystal field (CF) potential 
and a self-consistent octupolar interaction term 

XO(n(s)){0(n(s)))(T), (5) 

where n(s) represents one of the four inequivalent (1 1 1) 
directions, 0(n(s)) = J2 i=x y z Otm(s), A is a MF con- 
stant, and (0(n(s)))(T) is the self-consistent average 
value of 0(n(s)), which does not depend on s and is 
nonzero for T < To. The MF potential lowers the sym- 
metry of the Np Hamiltonian from Oh to D^d and splits 
the r 8 -quartet CF ground state into two singlets, T 5 
and Tq, and one doublet, T4, of D^- ^5 and Tq are 
complex-conjugate representations, which are degenerate 
for a time-reversal invariant Hamiltonian, whereas their 
degeneracy is lifted by the time-odd OP. By choosing the 
value of A which reproduces the observed To, the level 
sequence at T = is found to be singlet-doublet-singlet. 
The ground state has zero dipole moment and the sus- 
ceptibility saturates for T — > 0, as observed ||. These 
properties do not depend on the details of the CF Hamil- 
tonian (i.e. on the precise value of the CF parameter x 
or whether J-mixing is taken into account). 

This octupolar OP induces the observed triple-;? struc- 
ture of T 5 quadrupoles as secondary OP, with an ampli- 
tude cc p 2 in MF near To. The charge density is dis- 
torted along n(s), as quantified by the quadrupolar op- 
erator along n(s), $(ri(s)) oc a[3(n(s) ■ J) 2 — J(J + 1)] 
where a is a Stevens coefficient. $ is negative in the 
ordered phase, indicating a charge distribution "oblate" 
along n(s). This quadrupolar secondary OP is the quan- 
tity observed in the RXS experiment. 

The model we propose excludes a lattice distortion or 
a shift of the oxygen positions. The reduction of the local 
symmetry at the Np site leads to an electric field gradient 
along the (1 1 1) directions. This explains the line broad- 
ening observed below T in Mofibauer spectroscopy j|. 
The octupolar order breaks time reversal symmetry and 
thus allows the occurrence of interstitial magnetic fields 
as evidenced by ^SR 0. 

Finally, we remind the reader that we have inferred the 
triple-<f symmetry and octupolar order from the bound- 
ary conditions set out by previous experimental obser- 
vations, such as the absence of lattice distortions, the 
analogy to UO2, and from the excellent agreement of our 
data with the model. At the E2 edge, the principal OP, 
i.e. the T§ magnetic octupole, is expected to give a res- 
onant contribution to the scattered intensity jlTj], com- 
parable to that expected at the same energy from the 
secondary OP, the Tr, electric quadrupoles. However, no 
signal has been detected at the E2 edge, indicating that 
both contributions are below the sensitivity of the exper- 
imental device we used. At the El edge the intensity is 
almost entirely of quadrupolar origin, Indeed, the a — > a 



intensity has no magnetic contributions, and we can fit 
both the a —> a and the a — > n intensities using a single 
scaling factor, which should not be possible with a siz- 
able magnetic contribution in the a — > 7r channel. Such 
a contribution could originate from an ordered magnetic 
moment in the bulk, but Mofibauer spectroscopy estab- 
lish it unambiguously to be vanishing, or from dipolar 
order around defects, but this would give broad peaks 
which arc not observed. 

Whilst we have no direct evidence for a triple-<f mag- 
netic structure, and cannot directly prove the octupo- 
lar model, several predictions about the low temperature 
properties can be investigated. Octupolar order will split 
the Tg CF ground state into a singlet ground state and 
two excited levels. Inelastic neutron scattering H has 
shown an excitation near 6.5 meV, but the second one 
is expected to lie below 2 meV, outside the explored en- 
ergy range. Low temperature specific heat measurements 
should also show the characteristic signature of the CF 
spectrum. The integrated entropy from lowest tempera- 
tures up to the phase transition should be i?ln(4), with a 
significant contribution from a Schottky anomaly below 
8 K. If large single crystals would become available, oc- 
tupolar order could also be directly observed by neutron 
diffraction, with a form factor peaking at a finite value of 
the momentum transfer and reflecting the Fourier trans- 
form of the static magnetization density below To. 
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